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Abstract. We have theoretically investigated 3D focusing of a launched cloud 
of cold atoms using a pair of magnetic lens pulses (the alternate-gradient method). 
Individual lenses focus radially and defocus axially or vice- versa. The performance 
of the two possible pulse sequences are compared and found to be ideal for loading 
both 'pancake' and 'sausage' shaped magnetic/optical microtraps. It is shown that 
focusing aberrations are considerably smaller for double-impulse magnetic lenses 
compared to single-impulse magnetic lenses. An analysis of the clouds focused by 
double-impulse technique is presented. 



PACS numbers: 32.80.Pj, 42.50. Vk 

1. Introduction 

The field of atom optics pQ has undergone a dramatic expansion in the last two decades, 
largely as a consequence of the development of laser-cooling techniques 2 , and the 
routine production of atoms at microKelvin temperatures [3]. Paramagnetic cold 
atoms can be manipulated with the Stern-Gerlach force To date, the Stern-Gerlach 
force has been used to realise a variety of atomic mirrors for both cold [5] and Bosc 
condensed atoms This paper, however, concentrates on the formation of magnetic 
lenses for cold atoms. In comparison to a ballistically expanding unfocused cloud, a 
magnetically focused cloud can lead to density increases (or conversely temperature 
decreases) of many orders of magnitude. Applications include: atom lithography 7 ; 
transferring cold atoms from a Magneto Optical Trap (MOT) to a remote vacuum 
chamber of lower background pressure [8]; cold low-density atomic sources for fountain 
clocks 

The first demonstration of 3D focusing with pulsed magnetic lenses was conducted 
by Cornell et al. |1(J| . The group of Gorceix have made experimental and theoretical 
studies of cold atom imaging by means of pulsed magnetic fields ^^3- However, 
neither work addressed the optimum strategy for achieving a compact focused cloud, 
nor the limiting features for the quality of their atom-optical elements. 

Recently we provided a theoretical analysis of 3D focusing of weak-field-seeking 
cold atoms using a single magnetic pulse J3j- Lens designs for ID and 3D were 
presented that minimise aberrations due to the lens potential's departure from the 
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perfect parabolic case. Single-impulse 3D focusing has been experimentally achieved at 
Durham University, the results of which can be seen in a forthcoming publication ^3] . 

The scope of this paper is to investigate theoretically and numerically the limiting 
factors to the quality and size of the final image obtained in double-impulse magnetic 
focusing experiments; to identify the sources of aberration; and to discuss schemes 
for minimising their effect. We will show that both single- and double-impulse lenses 
yield a magnetically focused cloud with a bimodal distribution consisting of a highly 
diffuse outer cloud, as well as a core cloud which can be orders of magnitude denser 
than the initial atomic sample. This core cloud is therefore ideal for remotely loading 
tight traps with relatively small depth, e.g. miniature magnetic guides [1^], atom 
chips EI an d optical dipole traps JTj . 

The remainder of the paper is organised as follows: Section 2 outlines the theory of 
how to achieve the desired magnetic fields; Section 3 contains an analysis of magnetic 
imaging and minimising the final cloud size; Section 4 describes and contrasts the 
spatial performance of different magnetic lenses; Section 5 considers experimentally 
relevant parameters for alternate gradient lenses; Section 6 contains a discussion and 
concluding remarks. 

2. Alternate-gradient lens theory 

An atom in a magnetic field of magnitude B experiences a magnetic dipole interaction 
energy of U = —/i^B, where fi^ is the projection of the atom's magnetic moment onto 
the field direction. Provided that Majorana spin-flip transitions JS] occuring in field 
zeros are avoided and the rate of change of direction of the field is less than the Larmor 
frequency the magnetic moment adiabatically follows the field. 

The Stern-Gerlach force is Fsg = — V?7 = W(p(B). The ensemble (of alkali metal 
atoms) can be optically pumped into either a strong-field-seeking state with = /ib 
(where /ib is the Bohr magneton), or into a weak-field-seeking state with fj,^ = — /is- 
In low field, where the quantum numbers F and mp are good, these states are the 
stretched states \F = I + 1/2, mp = ±F). Atoms in these states have a magnetic 
moment which is independent of field, consequently the Stern-Gerlach force takes the 
simpler form Fsg = ±MbV_B — i.e. the focusing of the atoms is governed by the 
gradient of the magnetic field magnitude only. 

The choice of whether atoms in weak or strong-field seeking states are launched 
depends on the particular application. We discussed extensively in J2j the focusing 
of atoms in weak-field-seeking states. This was because single-impulse 3D imaging 
of strong-field-seeking states requires a maximum of the magnetic field in free space, 
which is forbidden by Earnshaw's theorem |19| . 

In this paper magnetic lenses centred on the point {0, 0, z c } are considered, with 
a second order magnetic field magnitude of the form: 



Bq and B2 are the bias field and the field curvature, respectively. Substituting this 
into the Stern-Gerlach force expression results in an atom of mass m experiencing a 
harmonic acceleration about {0, 0, z c } : 



twice the magnitude of, and opposite in sign to, the radial curvature, u> z 2 — — 2w r 2 . 




(1) 
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Note that lens curvature in all three spatial dimensions is reversed if the sign of either 
or B 2 is reversed. For simplicity from this point on, only the case /v = —/is i.e. 
weak- field-seeking atoms will be used and lens curvature is modified solely via B%. 

We refer to a lens with a field expansion of the form B as a one-dimensional lens, 
as it can be used either to focus axially (with simultaneous radial defocusing), when 
£?2 is positive, or to focus radially (with simultaneous axial defocusing) when B2 is 
negative. In order to achieve a 3D focus with the lenses of equation (j5J, an axially 
converging lens pulse must be followed by an appropriately timed axially diverging 
lens (or vice versa). This is referred to as the "alternate gradient" focusing method, 
and has the advantage of being able to focus both weak-field and strong-field seeking 
atoms. This method is used extensively in particle accelerators [2U, for focusing 
polar molecules an d is shown schematically in Figure ^ Useful parameters for 
describing the evolution of a Gaussian atomic cloud are the axial (c z ) and radial (ay) 
rms cloud radii (the standard deviations), as well as their aspect ratio £ = a* /ay. 



Figure 1. The principle of alternate gradient focusing. The upper image, 
Strategy AR, shows the evolution of the axial (a z , green) and radial (o> , black) 
cloud radii when an axially converging (radially diverging) lens precedes an axially 
diverging (radially converging) lens and leads to a sausage-shaped cloud (£ > 1). 
In the lower image, Strategy RA, the lens order is reversed, leading to a pancake- 
shaped cloud (£ < 1). The effects of gravity (in the —z direction along the coil 
axis) are not shown, but are included in simulations. Due to the time reversal 
symmetry of optics, the lens system also works backwards. 



As shown in reference |13| there exist optimal configurations for realising radial 
and axial focusing lenses. These are achieved with a pair of separated coaxial coils, 
where both coils have equal current with the same sense. An important lens parameter 
is the relative separation S of the coils in units of the coil radius. The harmonicity 
of a radially-focusing lens is optimized if S = 0.58 (red lenses in Figure ^); whereas 
the harmonicity of an axially-focusing lens is optimized if S = 2.63 (blue lenses in 
Figure |y. In the remainder of this work it is assumed that these optimized lenses are 
used. 




5=2.63 



Object (f=0) 
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3. Magnetic impulses and the ABCD formalism 

The separable equations of motion for an axially and cylindrically symmetric coil 
system lead to a lens that is harmonic in 3D with an acceleration given by equation 
J2J, which allows the motion in each cartesian dimension to be treated as a separate 
simple harmonic equation. The influence of a magnetic lens can be described by ABCD 
matrices, as outlined in |12l 113) . The initial and final position and velocity of an atom 
along a given Cartesian axis, say x, are related via the equation: 



Xf \ j A B \ ( x 



C D 



(3) 



(4) 



A 'thick' converging lens of strength uj (with Im(w) = 0) and physical duration r 
is actually equivalent to the 'thin' lens ABCD transformation: 

cos(wr) £sin(u;T) \ _ / 1 t'/2 \ / 1 \ f 1 r'/2 
-wsin(cjr) cos(wr) J V 1 / V C 1 / V 1 

where C{uj, t) = — uj sin(o;T), and is pre- and post-multiply by a translation matrix of 
half the effective pulse width t'(uj, t) = — tan(^). The notation of primes is used to 
denote times in the 'thin' lens representation. To model a diverging lens make the 
transformation uj — > ±\uj in equation I0J - i.e. C = ujsmh(ujT) and r' = ^tanh^p. 
The 'equivalent time' of the lens r' is not the same as the real experimental pulse 
duration of r. 



3.1. Double impulse magnetic lenses - the parabolic case 

A double lens system, see Figure El comprising lenses of strength and duration uji, t\ 
(starting after a time t\) and 0J2, T2 (starting a time £2 after the first lens) is modelled 
by using the following ABCD matrix sequence: 

U z)-(;?)U!)0*)U ?)(;?)• 

i.e. a t[ = ti + ir{ translation, then a strength C\ thin lens, a t' 2 — + t 2 + \t' 2 
translation, then a strength Ci thin lens followed by a t' 3 = T' — t[ — t' 2 translation, 
where Cj — —ojjSin(ujjTj) and rj = ^j-tan( ; ^ 2 -) for j E {1,2}. The total physical 
duration of the focusing, T, is fixed, and the effective total time of the double lens 
system is T = T - r x - r 2 + t[ + t' 2 . 

By multiplying the matrices of equation JSJ) together, the final ABCD system 
matrix is obtained. An image (i.e. a one-to-one map of position between the initial 
and final cloud) is formed if the condition B = is maintained. In this case the 
spatial magnification A is the inverse of the velocity magnification D\ a manifestation 
of Liouville's theorem. The cloud extent along x in a given plane is given by: 

al s ={Ao Xi f + {Ba Vi% )\ (6) 

where o~ Xi is the initial position standard deviation and a Vl . is the initial velocity 
standard deviation. An image is formed for the condition B = 0, but the smallest 
cloud size occurs when one minimises the product of the cloud extent for all 3 spatial 
dimensions (i.e. a rf 2 a Zf ). For single- and double-impulse lens systems, the cloud size 
at the image plane and the minimum cloud size do not correspond exactly, but they 
are usually very similar. In the rest of the paper we will consider the cloud size at the 
image plane (B = 0), and thus A corresponds to the magnification. 
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Figure 2. A schematic diagram showing the timing sequences for the two frames 
of reference, (a) shows the thick lens or lab frame with the time durations for 
each stage and (b) shows the mathematically equivalent thin lens representation 
that is used in the calculations. Note that the direction of time runs right to left 
so that it visually mirrors the system matrix layout in equation l5t . 



3.2. Solving the matrix equations 

The important entries of the system matrix in equation ^ are A and B: 

A=1 + (C 1 + C 2 ) (T' - t[) + C 2 (-1 + d (T' - t[)) t' 2 - C x C 2 t! 2 (7) 

B = T' + C 2 (T ; - t[ - t' 2 ) (ti + t' 2 ) + Citi (T' + C 2 (T' - t' 2 ) t' 2 - t[ (1 + C 2 t' 2 )) , (8) 

which are both second order in t[ and t' 2 (and hence also second order in t\ and t 2 ). 

To obtain an atom cloud which is focused in all 3 dimensions requires that 
the first lens is axially converging (radially diverging) and the second lens is axially 
diverging (radially converging), or vice versa. Moreover, the radial (subscript r) and 
axial (subscript z) spatial dimensions have different A and B coefficients. If the two 
axial lens strengths are u}\ z and uj 2z , then equation J5J) yields u>i r — iuji z /y/2 and 
uj 2r — iuj 2z /y/2. A 3D image is formed when equation JSJ) is set equal to zero for both 
the radial and axial directions. 

In reference |13| the density increase from a single-impulse isotropic 3D harmonic 
lens, A 3 /(I — A) 3 , was characterised by A, the equivalent time of the thin lens, t\\ 
relative to the total equivalent focus time T 1 . Note that for the anisotropic lenses 
in this paper the equivalent (i.e. thin lens) timing of a lens in the radial and axial 
direction is different. For this reason we will characterise alternate-gradient lensing 
with the parameters {Ai, A2} = {tx+Tx/2, ii+n +t 2 + r 2 /2}/T, corresponding to the 
mean times of the first and second magnetic impulses relative to the total experimental 
lensing time T. We use this labelling of {Ai, A2} if u>\ r is rea l (the first lens is radially 
converging), and we swap the definitions of Ai and A2 if io\ T is imaginary (the first 
lens is radially diverging). 

Modelling an experiment being conducted at Durham University |14| . we fix 
T = 212 ms. The cold atom cloud has an isotropic initial spatial and velocity 
distribution with ID standard deviations of o~r = 0.4 mm and ay = 4.4cm/s (i.e. 
a temperature of 20 /iK) respectively. The coils are assumed to have a 4 cm radius 
with 10,000 Amp-turn current in each coil. The two lens combinations in Figure^are 
shown in the table below with the resulting angular frequencies. 

For a range of values of T\ , and r 2 , we then solve the radial and axial simultaneous 
equations (jHJ (i.e. B r = 0, B z — 0) to determine t\ and t 2 . Although both B z and 
B r are quadratic in t\ and t 2 , substitution for either of these variables leads to a final 
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58irad s^ 1 


2.63 
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Table 1. The two different alternate-gradient strategies modelled. The uj's are 
the lens strengths, S's are the coil separations and £ = cr z /cr r is the cloud aspect 
ratio. 



sextic polynomial equation. This must therefore be solved numerically and leads to 
six solution pairs (ti,t 2 )- Only solution pairs with real times < < T satisfying 
the condition ti + t± + t 2 + t 2 < T are considered. The number of (ii,^) solution 
pairs as a function of T\ and r 2 is shown in Figure l^f a). These (iijta) solution pairs 
can then be used to calculate the relative increase in atomic density of a cold atom 
cloud. From equation @ the relative density increase of the image is thus: 

PaD = ^ _> 1 ( 9 ) 

((Aa fi ) 2 + (B r a v ) 2 )V(^R) 2 + (*W) 2 A r A z 

Where the arrow indicates the limit B = 0. The dimensionless relative density 
increases obtained for both strategies are shown in Figure |^b,c). These plots are 
then effectively combined in Figure G^d) by inverting t%, t 2 to find the relative density 
increase as a function of the parameters Ai, A2 which are the mean relative times of 
the radially diverging and converging lens. 



4. A measure of the quality of the focus 

The attributes of parabolic lenses are unimportant, unless it can be shown that 
experimentally realistic lenses are sufficiently parabolic for such an approximation 
to be appropriate. At the end of |13j there appeared to be a major difference between 
the parabolic approximation and real lenses. To some extent, this may have been due 
to the way the lens properties were measured. 

In numerical integration of the forces arising from the magnetic fields due to 
current loops (generated by the Biot-Savart law) was used, to track the trajectories 
of several (« 1000) simulated atoms. The initial positions and velocities of the atoms 
were randomly assigned, weighted according to isotropic Gaussian spatial and velocity 
distributions with ID standard deviations of an = 0.4 mm and cry = 4.4cm/s (as 
discussed in the previous section). The way in which the harmonicity of a lens was 
measured was to compare the expected harmonic focus size to the rms radii of the 
simulated atom cloud at the time of the harmonic focus. The important drawback of 
this rms approach is that the final location of atoms after a magnetic lens is highly 
nonlinear with respect to initial conditions. An atom with a velocity in the wings of 
the initial Gaussian distribution will experience highly anharmonic lensing, as it will 
be far from the centre of the lens during the magnetic impulse. Thus a few atoms can 
completely alter the rms width of the cloud. 

Another method to quantify the focus, adopted here, is to monitor the fraction 
of the atoms entering the focus region of a purely harmonic lens. The initial radial 
and axial cloud standard deviations are ctr, so the final standard deviations for a 
harmonic lens are ay — A r o~R and o~ z = A z o~r. By renormalising the dimensions 
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Figure 3. Image (a) shows the number of solution pairs (black=0, grey=l, 
white=2) for (ti,t2) as a function of t\ and T2 (t2 and n) in ms for Strategy 
AR (Strategy RA). The two dashed regions of the 'solution island' lead to the 
highest relative density increases, shown in (b) and (c). The relative density 
increase (equation J5J) if one images a cloud of atoms using: (b) Strategy AR, (c) 
Strategy RA. The maximum relative density increases are 1100 (320), for a £ = 17 
sausage (£ = 0.094 pancake) shaped cloud, for images (b) and (c) respectively. 
The results of (b) and (c) are combined in (d) , the relative density increase in terms 
of Ai and A2 (the mean times of the radially converging and radially diverging 
impulses relative to T). The points in images (b)-(d) are used later as a sample 
in simulations. 



so that the radial and axial dimensions are measured in terms of these final focus 
standard deviations, then a sphere with radius Rq = 1.53817 defined by 

L r^e 2 dr 1 , . 

J °° r 2 e ^dr ^ 

will contain half of the atoms of the focused Gaussian distribution. For numerical 
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simulations the fraction of atoms entering this harmonic focus is measured, and 
multiplied by twice the relative density increase of a purely harmonic lens (i.e. 
2A r ~ 2 A z ~ 1 ) to get a measure of the relative density increase afforded by a real lens. 

Note that the centre of the harmonic focus region, as well as the centre about 
which the rms radius is defined, is the final position of an atom initially in the 
centre of the Gaussian position and velocity distributions. This will lead to a slight 
underestimation in the density increase regardless of which way it is defined (it has 
been assumed the mean cloud position follows the initial centre-of-mass). In addition 
the density can increase if the best experimental 3D focus occurs at a time other than 
the best parabolic lens focus time, but the focus time has been allowed to vary in the 
simulations. For more details on the effects of gravity on launched atoms, namely that 
even the centre-of-mass atoms will experience time- varying radial and axial harmonic 
trap frequencies, see the Appendix. 

4-1- Single-impulse focusing revisited 

In light of the above discussion, 3D single-impulse focusing outlined in P2 is briefly 
revisited, in particular Strategy VI: the baseball coil system. This system consisted 
of a square baseball coil with side lengths w and a coil pair of radius w separated by 
w that were coaxial with the baseball coil axis. The current in the baseball coil was 
10,000 Amp-turns, and an isotropic lens is formed when the coil pair has a current 
of 1,541 Amp-turns. In this paper we have made the value w — 4 cm, leading to 
an angular frequency u> = 62rad/s in the harmonic lens. This is to provide a better 
comparison between single- and double-impulse techniques. 

In Figure El the relative density increase after the baseball lens is plotted in terms 
of the parameter A (the effective time of the baseball lens relative to T). The red 
dots correspond to the relative density increase using rms widths for the volume. The 
blue dots show the relative density increase as the fraction of atoms in the harmonic 
focus zone times the harmonic density increase. It is clear that we reach very different 
conclusions based on whether the rms radius of the focused atomic cloud, or the 
fraction of atoms which reach the harmonic focus are considered. 

By only looking at rms widths the optimum lens position occurs at A = 0.3, 
corresponding to a factor of 50 decrease in density. However considering the fraction 
of atoms in the harmonic focus zone times the harmonic density increase, the optimum 
position is now A = 0.9. The relative density increase is 2.3, which corresponds to 
0.3% of the cloud focused to a density 729 times greater than it was originally. 

The source of this discrepancy between methods can be seen in Figure 0] (b-c) 
where the non-Gaussian wings produce an over estimate of the rms cloud width. The y 
and z standard deviations for the Gaussian fits (black curves) are 730 /jm and 820 fim 
respectively; an almost isotropic distribution. 

4-2. Alternate gradient 

This section compares the alternate gradient numerical simulations with the purely 
harmonic lenses of section[3] The (ti, t%) sample co-ordinates illustrated in FigureOJb- 
c) are used, in order to run numerical simulations for the relative density increases 
illustrated in Figure We have not (cf. Figure 0J used the rms volume of the cloud 
to show the relative density increases as these result in extremely low relative density 
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Figure 4. Image (a) shows the relative density increase (with error bars, for 
a 5000 atom simulation) based on: a purely harmonic lens (black curve), the 
fraction of atoms in a real lens arriving at the harmonic focus region (blue dots), 
the ratio of rms cloud volume before and after a real lens (red dots). Images (b) 
and (c) show the spatial probability distributions at the focus (black dots) in the 
y and z directions, respectively, where the non-Gaussian wings of the distribution 
can be clearly seen. These two distributions are taken from the A = 0.5 lens 
simulation, in which case 16% of the atoms are in the harmonic focus region. The 
Gaussian fits in (b) and (c) (as well as the x distribution) have an area of PS 70%. 



increases (typically 1CP 5 in (a) and 10~ 3 in (b)) that would reduce the contrast in 
Figure EJ 

For both strategics the numerical simulations trace the shape of the analytical 
relative density increase although aberrations result in reduced increases. The 
maximum relative density increases in (a) and (b) are 186 and 50 respectively. This is 
a marked improvement on the single-impulse focusing, however the cloud distribution 
is no longer isotropic. The harmonic focus aspect ratio has a range 12 < £ < 17 in (a) 
and 0.074 < £ < 0.095 in (b). 

In certain applications, for example microtrap loading and lithography, the 
sausage-shaped distribution with its reduced radial spread could be beneficial. Figure|S] 
(c) plots the distribution of a cloud focused via Strategy AR. The standard deviations 
for the Gaussian fits to the core of the x, y, and z distributions are 56 /mi, 56 fim 
and 850 [im; an order of magnitude reduction in the radial direction compared with 
single-impulse focusing. 

Given comparable lens dimensions and strengths, double-impulse magnetic 
focusing is far superior to single-impulse magnetic focusing in terms of the relative 
density increases that can be achieved by a fraction of the atoms. This result is in 
stark contrast to the relative rms density increase of the entire cloud, which would lead 
to the opposite conclusion. This spatial focusing would find applications in lithography 
or sending the atomic cloud through micro-sized apertures. 

5. Physical properties of clouds focused via the alternate gradient method 

For many experiments it is also important to consider changes to the velocity 
distribution and hence the collision rate and phase-space density. Unless the collision 
rate of a gas is sufficiently high, then Bosc-Einstein condensation via sustainable 
evaporative cooling is impossible. Changes in the collision rate and phase-space 
density of the focused cloud are complicated by the fact that alternate-gradient lensing 
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Figure 5. Images (a) and (b) use the (t\,T2) co-ordinates illustrated in 
Figure |3b) and (c) to show the relative density increase for alternate gradient 
lensing Strategies AR and RA, respectively. There were 1000 atoms in the 
simulation and relative density increases are shown for a pure harmonic lens (black 
curve), as well as the relative density increase for the fraction of atoms in a real 
lens arriving at the harmonic focus region (blue dots with error bars). In image 
(c) the strong spatial bimodal nature of the x focus for the leftmost point in (a) 
is clearly seen on a log scale. The Gaussian fit (with a x = 56 fim) contains 49% 
of the 3000 simulated atoms used in (c). 



automatically leads to an anisotropic focused distribution both in space and in velocity. 

At the B T>Z = focus of the cloud the ABCV ryZ matrix has radial and axial spatial 
widths a T — (r 2 ) 1 / 2 = A r OR and a z — (z 2 ) 1 / 2 = A z or respectively. The radial and 
axial velocity widths are given by: 

\]{v r ,z 2 ) = y/k B T riZ /m = ^C r /<7 R 2 + a v yA r , z 2 , (11) 

where T TtZ is the atomic cloud temperature. In the limit C TtZ a R <C o~v/A r , z (which is 
not always the case) , the velocity width of the focused cloud is inversely proportional 
to its spatial width and \J (v r , z 2 ) — &v /A r , z . In Figure El phase-space plots of the AR 
and RA strategies from the rightmost points of Figures |S] (a) and (b) were generated 
with a 30,000 atom simulation. The effect of aberration is clearly seen when comparing 
the purely harmonic lenses (subscript H) and the Monte Carlo simulation with the 
full magnetic fields from realistic coils (subscript MC). The plots also demonstrate the 
inversely proportional relationship between spatial and velocity widths. Furthermore, 
in the Monte Carlo simulations there is much stronger correlation between position 
and velocity. 

The anisotropic temperature in the focused cloud means that captured atoms in 
a cylindrically symmetric harmonic trap with radial and axial frequencies v r and 
v z , will rethermalise to an isotropic temperature via elastic collisions. The total 
(potential+kinetic) energy of the focused cloud is equated with that of a 3D harmonic 
oscillator at equilibrium temperature T (i.e. hhsT average energy per atom for each 
space and velocity dimension). The equilibrium temperature of the focused cloud after 
thcrmalisation is thus: 

T = (4^ 2 (2^V) + ^ V» + (2<tv 2 ) + (O)) • (12) 

In many experiments other physical properties of the atomic cloud are of interest: 
the atomic density n cx l/(<7 r 2 a z ), collision rate 7 cx nVT and phase-space density 
PSD cx n/(T r T z 1/2 ). In order to minimise the loss of phase-space density during 
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Figure 6. Phase-space plots of the AR and RA strategies from the rightmost 
points of Figures [F] (a) and (b) respectively. Both a harmonic lens calculation 
(subscript H) and a Monte Carlo simulation with the full magnetic fields from 
realistic coils (subscript MC) are plotted. In the ARmc pl°t there are 45% (x — v x ) 
and 42% (z — v z ) of the initial 30,000 atoms present. In the RAmc plot there are 
61% (x — v x ) and 24% (z — v z ) of the atoms present. The effects of aberration 
are clearly seen when comparing the harmonic and realistic magnetic coils. 



thermalisation of the focused cloud, one can show that it is best to choose v r and v z 
such that the potential energy is equal in all spatial dimensions and the total potential 
energy is equal to the total kinetic energy of the cloud. 

The physical properties of the focused atomic clouds of Figure are displayed 
in Table |21 The relative values for the density n, collision rate 7 and phase-space 
density can be converted into absolute values using typical initial experimental values 
(n) = 10 10 cm~ 3 (an atom number N = 3 x 10 7 in the unfocused cloud), 7 = 1 Hz (for 
87 Rb with the low-temperature collision cross-section a = 8ira 2 = 8x 10~ 12 cm 2 ) and 
phase-space density PSD=2 x 10~ 6 . Aberrations will effectively be smaller if a large 
MOT is used with the same density (e.g. N = 10 9 , ctr = 1.3 mm) and temperature as 
Table |21 and a trap with weaker frequencies to catch the atoms. 

Note that the Monte Carlo relative density increases in Table [2] are lower than 
those in Figure [5] by a factor of « 2, but the fraction of atoms in the focus are higher 
by a factor of f» 2. In Figure [S] the density is estimated by measuring the fraction of 
atoms arriving at the harmonic focus. Here, a 6-dimcnsional Gaussian phase-space fit 
to the narrow central peak of the bimodal focus was made to explicitly obtain ay, o~ Zl 
T r and T z . The fraction of atoms in this Gaussian focus was used for P here. 

An interesting result of Tablets that the aberrations of 'real' lenses work to our 
advantage, to some extent, in that the atoms can be loaded into a trap with a shallower 
depth than atoms focused by a purely harmonic lens, and phase-space density loss is 
reduced during rethermalisation. This is due to the reduced anisotropy of the spatial 
and velocity distributions at the focus. The focused atoms have a relatively high 
temperature and one needs a trap depth of ss 10 mK to trap the focused atoms. For 
an atom with a magnetic moment of one Bohr magneton /zg, this corresponds to a 



Double-impulse magnetic focusing of launched cold atoms 



12 





t 


P 


oy iZ (/im) 


T r . 2 ( M K) 


Vr,z (Hz) 


n 


7 


PSD 







1 


400, 400 


20, 20 


17.4, 17-4 


1 


1 


1 


AR H 


T 


1 


15.2, 257 


13900, 60.9 


12000, 47.3 


1080 


23200 


0.892 




Tr 


1 


15.2, 257 


9260, 9260 


9850, 584 


1080 


23200 


0.108 







1.000 


400, 400 


20, 20 


17.4, 17.4 


1 


1 


1 


ARmc 


T 


0.308 


23, 500 


2720, 57 


3530, 23.5 


74.5 


713 


0.325 




Tr 


0.308 


23, 500 


1830, 1830 


2900, 133 


74.5 


713 


0.085 







1 


400, 400 


20, 20 


17.4, 17-4 


1 


1 


1 


RA h 


T 


1 


128, 12.2 


197, 21600 


170, 18700 


319 


6110 


0.988 




Tr 


1 


128, 12.2 


7320, 7320 


1040, 10900 


319 


6110 


0.046 







1.000 


400, 400 


20, 20 


17.4, 17-4 


1 


1 


1 


RAmc 


T 


0.305 


200, 20 


105, 1670 


79.7, 3180 


24.4 


137 


0.509 




Tr 


0.305 


200, 20 


627, 627 


195, 1950 


24.4 


137 


0.139 



Table 2. Physical properties of the two different alternate-gradient strategics 
modelled, AR and RA from the rightmost points in Figs.|F](a) and (b) respectively. 
Subscripts H and MC respectively denote a simulation with purely harmonic lenses 
and a Monte Carlo simulation with the full magnetic fields from realistic lens coils. 
The measured parameters are: fraction of atoms in the Gaussian focus P (see 
text), radial/axial cloud radius <T r ,z, radial/axial temperature T r ,z, radial/axial 
trap frequency u T}Z , relative density n, relative collision rate 7 and relative phase 
space density PSD. The effective trap frequencies for the initial (t = 0) and focused 
(t = T) cloud (italicised) are equilibrium values based onv = y T/m. The 
actual frequencies of the trap the atoms are loaded into at t = T are denoted in 
the t = Tr lines (the cloud properties after T plus the thermalisation time). 



magnetic trap depth 150 G. 

6. Discussion and conclusion 

The main application of interest for the magnetically imaged atoms will be loading 
a magnetic microtrap or optical dipole trap, for which alternate gradient imaging is 
well-suited. If the trap that is being loaded is harmonic, with a large capture volume, 
then the rms size of the cloud will be linked to the equilibrium temperature after clastic 
collisions rethermalise the initially bimodal image distribution. In order to keep the 
high density core of the atomic cloud, the high energy atoms must be removed on a 
time scale that is rapid compared to rethermalisation - this could be achieved with 
strong RF evaporative cooling or by shining resonant dark SPOT beams [22] at the 
focal region. A trap with a small capture volume, e.g. an atom chip |16j or a focused 
optical dipole beam trap JTj, is ideal as only the high density core of atoms will be 
captured in the trap. 

The timescale needed to remove the shell of hot diffuse atoms from the central core 
will actually be much longer than the time it takes for the dense core atoms (which are 
anisotropically focused in space and velocity) to thermalise to a uniform temperature 
in 3D. The higher the anisotropy of the focused atomic cloud, the more phase space 
density is lost during the thermalisation. This loss could be mitigated to some extent 
if the radial and axial trap frequencies were varied to remain in equilibrium with the 
changing radial and axial cloud temperatures during the (short) thermalisation time. 

In this paper we have used experimentally realistic parameters to compare the 
limiting focal size of a launched cold cloud of weak-held seeking atoms subject to 
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either a single or double magnetic lens impulse. The ABCT> matrix formalism was 
convenient for giving an estimate as to the parameters needed for magnetic focusing, 
but numerical simulations were necessary to detect the effects of aberrations in real 
magnetic lenses. If one wishes to minimise the rms image volume of a launched cloud 
then a single-impulse lens is preferable. If, however, one can selectively capture the 
central core of the bi-modal image, a double-impulse (alternate gradient) lens can lead 
to orders of magnitude relative density increases for both pancake- and sausage-shaped 
image clouds. Although we have only considered cold thermal atomic clouds in this 
paper, the effects of aberrations will also be important for tight focusing of coherent 
matter waves EH] ■ 
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Appendix: modelling time varying lens strengths 

This paper has only discussed parabolic lenses with constant strength to, pulsed on 
with a top-hat pulse of duration r. In experiments the lens strength is a function 
of time, partly because the current in a real coil is not a top-hat pulse, and partly 
because the centre-of-mass of a launched atomic cloud changes as it goes through a 
lens and thus (to second order in position) will experience a time-varying parabolic 
lens. In practice it was unnecessary to adjust the timing or lens coil positions to allow 
for these effects in the simulations, however we briefly discuss ways around this issue 
should it become problematic. 

A harmonic lens with arbitrary time variation uj(t) (where uj(t) and its derivatives 
are zero outside the experimental pulse time t\ < t < ti+Ti), is equivalent to an infinite 
product of infinitesimal 2x2 translation and thin lens matrices, resulting in a single 
2x2 matrix that is itself independent of the initial velocity and position of an atom. 
On solving y"[t] = — uj 2 [t]y[t] from t = t\ to t = ti + ti, with the initial conditions 
{y[ti],y'[ii]} = {dy,0} to get {y[t 1 +n],y'[ti + n]} = {A,C}dy, then use initial 
conditions {y[ti}, y'[ti]} = {0,dv} to get {y[t\ + Tx],y'[ti +ti]} = {B,V}dv, (with 
small values for dy and dv) results in the general ABCT> matrix for any initial position 
and velocity from (numerically or analytically) solving the differential equation for 
only two different initial conditions. As this ABCT> matrix has determinant 1 (it is a 
product of determinant 1 matrices) it can be expressed as a (translation matrix)- (thin 
lens)- (translation matrix) combination. 

One can then use an iterative 4D Newton-Raphson method with four input 
parameters (the z positions of the two alternate gradient lenses and the times t\ 
and £2) such that the z centre-of-mass velocity of the atoms is not altered by either 
lens, and B r and B z are identical to zero. 
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